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In addition to the scalar Shakura-Sunyaev a ss turbulent viscosity transport term used 
in simple analytic accretion disc modeling, a pseudoscalar transport term also arises. The 
essence of this term can be captured even in simple models for which vertical averaging is 
interpreted as integration over a half-thickness and one separately studies each hemisphere. 
The additional term highlights a complementarity between mean field magnetic dynamo 
theory and accretion disc theory treated as a mean field theory. Such pseudoscalar terms have 
been studied, and can lead to large scale magnetic field and vorticity growth. Here it is shown 
that vorticity can grow even in the simplest azimuthal and half-height integrated disc model, 
for which mean quantities depend only on radius. The simplest vorticity growth solutions 
seem to have scales and vortex survival times consistent those required for facilitating planet 
formation. Also it is shown that when the magnetic back-reaction is included to lowest order, 
the pseudoscalar driving the magnetic field growth and that driving the vorticity growth will 
behave differently with respect to shearing and non-shearing flows: the former can reverse 
sign in the two cases, while the latter will have the same sign. 
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1. Introduction 

Analytic accretion disc theory and mean field magnetic dynamo theory are often thought 
of disjointly, but the two are intimately related in certain important senses and exploring this 
relation has some interesting consequences. In traditional a analytic accretion disc theory, 
(e.g. Shakura & Sunyaev 1973; Pringle 1981; Balbus-Hawley 1998) the Reynolds+Maxwell 
stress tensor representing a correlation of turbulent velocities, is replaced by a ss c^, where c s is 
the sound speed. In the Navier-Stokes equation, this amounts to replacing the microphysical 
viscosity by a turbulent viscosity, v = a ss c s H, where H is the scale height. The flow variables 
then being solved for must be interpreted as mean fields, with gradient scales larger than 
those of the turbulent motions providing the viscosity. (Consequences of fluctuations for 
variability were studied in Balbus et al. 1994; Blackman 1998; 2000a). The continuity, 
radial momentum, angular momentum, and energy equations are solved, often with the 
assumption of axisymmetry, for radial solutions. The magnetic field is often added in as a 
pressure, rather than coupling in the magnetic induction equations (e.g. Narayan & Yi 1994; 
Narayan et al. 1998). 

Indeed it is widely known that the a viscosity approach is incomplete. For more realistic 
discs, the dynamical evolution of the magnetic field must be coupled to the fluid equations. 
The problem is fully non-linear, and ultimately requires non-linear MHD turbulence simu- 
lations understanding magneto-shearing instabilities etc. (c.f. Balbus & Hawley 1991,1998, 
Brandenburg et al. 1995; Miller & Stone 2000). Nevertheless, the a viscosity approach 
provides a useful phenomenological framework for exploring solutions. This is where the 
relationship to mean field dynamo theory enters. 

As traditional accretion disc theory ignores the induction equation, traditional kinematic 
mean field dynamo theory solves only the magnetic induction equation (e.g. Moffatt 1978; 
Parker 1979; Krause & Radler 1980; Radler 1999) taking the velocity that appears in that 
equation as given and not subject to the influence of the field. Thus traditional accretion 
disc theory and kinematic magnetic dynamo theory are complements of each other. Mean 
field dynamo theory in its analytic regime, kinematic or dynamic, is also subject to the 
complementary criticisms of accretion disc theory in that rigorously the fully dynamical set 
of equations have to be solved and the physics of MHD turbulence must be understood. For 
a sheared accretion disc, the problem of accretion and the problem of field growth are one 
fully coupled problem. Disc simulations with appropriate boundary conditions may show 
some evidence for mean field dynamo operation in a turbulent accretion disc (Brandenburg 
& Donner 1997), though higher magnetic Reynolds number simulations are needed. Note 
also that the mean field dynamo is not to be confused with turbulent amplification (= small 
scale dynamo). The former models an inverse cascade of magnetic helicity (e.g. Pouquet et 
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al. 1976), not simply field line stretching on the scale of the turbulent velocities. In addition, 
the mean field dynamo is often invoked to generate a mean magnetic flux, not just a mean 
field energy. 

The formalism of mean field theory as applied to the induction equation leads to both a 
viscosity term that corresponds to the Shakura-Sunyaev a viscosity term, and an additional 
pseudoscalar term which can be similarly parameterized. As applied to the Navier-Stokes 
equation, the same formalism leads again to a Shakura-Sunyaev a ss type viscosity term, 
and also to a pseudoscalar helicity term. The latter would vanish if the vertical averaging is 
interpreted as integration over the entire disc height, but does not vanish when the integration 
is taken over a half scale height (or less) and one considers each hemisphere independently. 
The pseudoscalar term has been recognized elsewhere (e.g. Moiseev et al 1983; Frisch et al. 
1987; Khomenko et al 1991; Kitchatinov 1994ab; Tanga et al. 1996; Blackman & Chou 1997). 
Here it will be shown that it can be treated on par with the Shakura-Sunyaev viscosity. The 
pseudoscalar term can lead to vorticity growth in discs. 

Enstrophy exhibits an inverse cascade in 2-D turbulence (Kraichnan & Montgomery 
1980). The growth of vorticity in primarily 2-D rotating fluids has been seen in nature (e.g. 
Jupiter, c.f. Ingersoll 1990; Marcus 1993) as well as in simulation (Marcus 1990; McWilliams 
1990) and experiment (e.g. Sommeria et al. 1988). Statistical mechanics approaches have 
been successful in modeling this (e.g. Chavanis & Sommeria 1998). Vorticity growth in 
sheared accretion discs is less well studied (reference added in proof: Godon & Livio 2000), 
even though vortex evolution has been studied (e.g. Adams & Watkins 1995; Godon & Livio 
1999; Chavanis 2000). It will be later shown that an instability results in a sheared disc even 
when mean quantities are integrated to be only a function of radius in each hemisphere. 
Note here that the use of a pseudoscalar requires the turbulence to be 3-D even if the mean 
fields are taken to be 2-D. 

In sections 2 and 3 the mean vorticity and mean magnetic field equations are derived with 
the back-reaction between the two included to lowest order, under the assumption that 1st 
order cross-correlations between fluctuating fields and velocities are small. In section 4, the 
resulting diffusion and pseudoscalar transport coefficients a parameterized to be proportional 
to the sound speed and, growing solutions for vorticity are derived. In section 5 the relation to 
vorticity growth for planet formation is discussed. Why the pseudoscalar growth coefficient 
for vorticity can have the opposite sign to that driving mean magnetic field growth is also 
discussed. Section 6 is the conclusion. 
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2. Derivation of Vorticity Equation 

The Navier-Stokes equation with B-fields is given by 
dtv = v x (V x v) - Vp/p - Vv 2 /2 + z/V 2 v + (C + |^)V(V • v) + m + F(x, t) + V0, 

(1) 

where v and £ are the constant viscosity and second viscosity and p is the pressure and 
m = (B • VB — |VB 2 )/47rp. The V0 includes the potential force of gravity. The vorticity 
equation, 

d t u> = V x (v x u) + vV 2 u + Vp x Vp/p 2 + V x F(x, t) + V x m, (2) 

where u> = V x v, is obtained by taking the curl of (|l|). Hereafter we ignore the microphysical 
viscosities. 

Now divide quantities such as u>, v, etc. into mean (indicated by an overbar or ( )) and 
fluctuating (indicated by prime) components v = v + v' and w = u> + u/ and respectively. The 
interest is in the application to a thin disc. I take the average to mean full integration over 
azimuth, vertical integration over the top half of the disc only, and local radial averaging, 
leaving mean quantities only a function of radius in each hemisphere, considering the top and 
bottom hemispheres of the disc separately. This is important because pseudoscalar averages 
would vanish if the average is taken over the whole disc. 

The presence of F(x, t) represents any function (e.g. forcing) not included in the other 
terms. It will drop out later. Note that ultimately, pseudoscalar generation results from 
a vertical density gradient combined with the shear and underlying rotation. Later I will 
parameterize the pseudoscalar in a similar way that Shakura & Sunyaev (1973) parameterize 
the viscosity. 

It is assumed that derivatives with respect to x or t obey dt tX (XiXj) = (dt, x (XiXj)) and 
(XiX'j) = (Reynolds relations ( |Radler 198(J| )), where Xj = Xi+X[ are components of vector 
functions of x and t. For the spatial mean, defined by pQ(x, £)) = |C|~ 3 J^-l ^( s ? t)d 3 s, the 
relations hold when the averaging is taken over a large enough scale, such that I <C \(\ < L, 
where L ~ v/Vv, and i ~ v'/Vv'. 

Subtracting the mean of (]I|) from itself, and assuming V0' = 0, gives 
d t V = (V ■ Vv') - v' • Vv' - v • Vv' -v'-Vv 

-Vp'/p - V(v' ■ v) + F'(x, t) + m' ! ° ' 

ignoring the viscosities. To proceed, (Q) needs to be averaged. To simplify, assume that 
the density has only a mean time independent spatially varying quantity, and ignore its 
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fluctuating gradients. Then, since the means are only functions of the radial coordinate, the 
pressure term drops out of the vorticity. (The mean pressure thus satisfies the barotropic 
equation of state p = p(p).). This gives 

d t u3 = V x (v'xu'} + Vxm + w-Vv-wV'V, (4) 

where I neglected terms second order in time-varying mean quantities and and d t indicates 
working in the Local Standard of Rest (LSR) frame, i.e. that which moves with V. I have 
also assumed F = F'. Subtracting (Q) from (0) (ignoring the pressure term for reasons 
described above) gives 

d t J = J • Vv - JV • v + zu • Vv' -wV-v' 

-v' • VaJ + J ■ Vv' - u/V • v' - v' ■ W - V x (v' x J) + V x m' (5) 
+V x F'(x,t). 

For u; to grow, the Vx terms in (|J) must be non- vanishing. 

Following previous work (c.f. Blackman & Chou 1997; Field et al. 1999), I expand the 
turbulent quantities on the right of Eqs. (|J) to linear order in Vv using the equations for 
the fluctuating fields here it is assumed that correlations of zeroth order quantities can be 
functions of radius. To find the lowest order terms, assume weakly anisotropic turbulence: 
terms linear in the mean shear contribute, but their averaged th order coefficients, are taken 
to be isotropic. These coefficients can be reflection asymmetric and radially dependent. 

In order to "ignore" the terms which are products of one zeroth order turbulent quantity 
with one 1st order turbulent quantity (as I will do here) one assumes that these terms are 
small compared to the associated terms that involve products of one zeroth order quantity 
with one mean quantity. For example, one would assume that the 1st term on the right of 
(H) is greater than the sixth term on the right of @. This weaker than the "usual" first 
order smoothing approximation in that the present requirement amounts to a/^ < Vv and 
v^ < v, rather than the usual J < Vv or v < v. The latter two conditions are stricter as 
they apply to the total fluctuating quantities rather than only the anisotropic part. 

An alternative potential justification for dropping the troubling terms that are products 
of one zeroth order turbulent quantity with one 1st order turbulent quantity, is that they 
themselves likely have positive and negative contributions which might nearly cancel: In a 
steady state, the turbulent quantities are balanced by input and cascade. This motivates a 
possible replacement of e.g. the sum of 6th 7th 8th and 9th terms by (x — £) w > where x an d 
£ are positive, and their difference represents the combination of growth and decay. (The 
9th term is actually irrelevant as it vanishes when correlated with a fluctuating quantity, 
which is the only context in which it will enters.) Ignoring these terms would then amount 
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to the assumption that ( — £ is small. Such an approach avoids comparing each individual 
"offending" term to those dependent on the mean fields, since here the offending terms would 
cancel themselves. 

Working in the LSR frame is also important to emphasize here. Although the turbu- 
lence in a sheared disc is more than weakly anisotropic, by working in the LSR frame, the 
anisotropy then manifests through terms like the first in ([5]) oj'^ ■ Vv rather than through 
terms like to v • Vw'^. To assess the implications, terms like v'(°) • Vv must be compared 
to terms like v'(°) • Vv'' '. For accretion discs whose turbulence is ultimately driven by a 
magneto-shearing instability, the ratio of these terms is of order * a But /It ~ f2 for 

(v^> /It) 

the magneto-shearing instability, so the anisotropy is of order 1 in the LSR frame, compared 
to >> 1 in the lab frame. Indeed, even in the LSR, this still means that the anisotropy 
should be considered to more than linear order, but because it is "only" of order 1 we expect 
qualitative similarities of the results to the expansion to all orders. Blackman (2000b) gave 
a restricted approach to treating this shear anisotropy to all orders. 

Using the formal solutions for the turbulent field u/(t) = <J{t = 0) + / dfu'dt', and using 
times appropriately chosen such that the correlation (v'(i) x oj'(0)) ~ 0, we get 

(V x u/) (1) = (v'(°)(t) x / ' d t ,J^dt') + (J* dvv'Wdt' x u/(°)(t)>, (6) 

The calculation of these averages requires Eqs. (^) and (||) for the time integrands, invoking 
the approximations discussed above. Using (^) also requires an expression for the pressure, 
which arises in (^) via the term 

7TV (0) (*) x tW^dt'), (7) 
J o 

where the p is pulled out under the assumption that it is time independent. Using isotropy, 
homogeneity, and Reynolds rules, Blackman and Chou (1997) showed that terms of the form 
(0) vanish in the derivation of the mean field equations to the first order considered. However 
in the present case, the averages of statistical correlations can be functions of radius. Thus 
several terms dropped from Blackman & Chou (1997) must be considered. Following that 
approach, the energy equation can be used to obtain 

V/ 1} = V/ 1} (0) - V / (v ,(0) -Vp + v-Vp ,{0) (8) 
+7p'(°)V- v + 7 pV- v /(0) )dt', 

where 7 is the adiabatic index and we must now relax the V-v"(°) constraint. To simplify the 
analysis, I have ignored any radiative contributions to higher than zeroth order quantities. 
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The pressure dependent contribution to (v' x lv')^ can then be written 




= 2 -f J*'dt'[d k p((Vv') 2 )^ - d&idjvfttyW - TD k (W (0) ■ V x v'(°)>], 



(9) 



where it is assumed that mean fields vary on time scales longer than the fluctuating fields, and 
Vp'^(t = 0) is uncorrelated with or w'(t). The vanishing of terms to get to the last equality 
follows from careful application of isotropy (i.e. rank 2 and rank 3 averaged tensors of 
fluctuating th order quantities are proportional to and respectively), but homogeneity 
of the th order turbulence has not been used (i.e. di(XjXk)^ = 0) nor has V ■ v' = been 
used. Had homogeneity been used, then the first term on the right of the last equality would 
not have survived. Note however, that when put inside the curl of (B), the first two terms 
on the right of (^) vanish. This is because the curl of the pressure gradient vanishes, and all 
gradients are in the radial direction, so the cross product of gradients vanishes. Finally, the 
last term in @ vanishes from use of Reynolds rules and isotropy, which is seen from using 
the chain rule with the gradient on p'^ . The pressure does not seem to contribute to lowest 
order under the approximations used herein. 

Collecting all of the above, crudely approximating time integrals by factors of the cor- 
relation time r c (Ruzmaikin et al. 1988), and freely employing Reynolds rules and incom- 
pressibility of fluctuating components, 



where it is assumed F = F^ ), and the microphysical viscosity is ignored. The coefficients 
are 



V x (v' x u/) (1) = V 



x a u> — V x (PqV x 13 



(10) 



where b'^ = B'(°> /47rp 




(11) 



d t (jj 



V x a u — V x (/3 Q V x w) + uj ■ Vv — wV • v. 



(12) 



This equation presumes that 1st order cross correlation terms vanish, that is 



(b'(°)-v'(°)) = (a;' (0) -b'( )) = (u;' {0) -Vx 



b'<°>> = 0. 



(13) 



When this assumption is not made, then the mean magnetic field and mean vorticity field 
equations are coupled (Blackman & Chou 1997). Note that since the averaging is such that 
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the mean fields are only functions of radius, the penultimate term in ( |12D vanishes. The last 
term vanishes by using the continuity equation for a steady mean density (which has only a 
mean component) which gives pV ■ v = — v • Vp, and then recalling that we are working in 
the LSR frame. Thus 

d t Z3 = V x a uj - V x (/3b V x w) (14) 

is the equation to be solved. Before proceeding to do so, I derive the magnetic field growth 
equation for comparison. 



3. Derivation of Magnetic Field Equation 

As mentioned above when cross correlations between velocity and magnetic field com- 
ponents are ignored, the mean field equations for the magnetic field and the vorticity field 
decouple for the case in which the gradient of the fluctuating components of the density are 
ignored. The induction equation for B is then 

d t B = V x (v x B) + z/ M V 2 B, (15) 

Similarly, the equation for the mean B-field, derived by averaging ([15]) is 

d t B = V x (V x B') + B Vv, (16) 

where the resistivity has been ignored and V ■ v has been assumed. Subtracting (|I6"D from 
(|T5|) yields the equation for the fluctuating B-field 

d t B' = B' • Vv - v VB' + B • Vv' — BV • v' — v' • VB+ 

B' • Vv' - B'V • v' - v' • VB' - V x (v' x B') ^ ' 

again ignoring the resistivity, and assuming V ■ v = 0. Following the same procedure to first 
order in B and v for (v' x B') in ( |i~6|) that was followed above for (u/ x v') assuming (|13|) 
holds, the mean field induction equation becomes 

d t B = V x a m0 B - V x (PrM V x B) + B • Vv, (18) 

where the coefficients are 



a m0 = (r c /3)(2(B'(°) ■ V x B'(°))/47rp - (v'W • V x v'^)) 
Ano = (r c /3)(2(B'(°) • B'(°))/4vrp + (v'^ ■ v'^)). 



(19) 



and we have used p = p. 
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4. Instability growth 

Since the mean quantities are taken to be only a function of radius in each hemisphere, 
it is convenient to break up the vorticity equation into toroidal and poloidal components. 
Recalling that we are working in the LSR frame so that we can ignore terms that include 
factors of v without derivatives, the two resulting equations that we need to solve are then 
of the form 

d t uj^ = «oV xu P + Va x up + /3 V 2 o^ - V/3 x (V x ufy) (20) 

and 

dtv<t> = a W^ + /^oV 2 ^, (21) 

where the subscript indicates the toroidal component and P the poloidal component. 

Writing these equations for their components in cylindrical coordinates, and recalling 
that only the radial spatial derivatives of mean quantities contribute, gives 



^ - ^d r a - (d r a)(d r vt) + ^(d r (3 ) + (d r fa)(dj3+) 



(22) 



and _ 

d t v^ = a Q W^ + PodrVc/, + —d r v^ (23) 



We look for solutions to fl22|) and (^) of the form 

= W^ (r)e ifc ^ + "' and ^ = %)(r)e^ +nt , (24) 

where k r , is the radial wave number, Uk is the Keplerian speed, v^o(r) = Vk(r) oc r _1//2 , and 
w^o (r) oc r~ p . We can take p = 3/2 since all velocities will scale with the Keplerian speed, 
and the vorticity scales with the associated curl. 

Since the correlation coefficients can also depend on r, this dependence must be ad- 
dressed. From (??) it is evident that (3q is the turbulent viscosity which can be parameterized 
in the Shakura-Sunyaev (Shakura & Sunyaev 1973) form, that is 

Po = a ss c s h, (25) 

where c s is the sound speed, h is the disc half thickness, and a ss is the Shakura-Sunyaev 
parameter. Similarly, «o can be parameterized as 

«o = qoissCs = qflo/h, (26) 
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where q satisfies — 1 < q < in the top half of the disc, and < q < 1 in the bottom half. 
Since c s ~ hvk/r, for h oc r, (3 oc r -1 / 2 and «o oc r~ 3 / 2 . Taking into account the radial 
dependences, and plugging ( ]2"3| ) into ( |2"2"D and ( |2"B| ) gives 



2 3/3 5 ifc r /3 

n + /3 k r - -— + 

2r 2 r 



ou - 



+ 7:^7 + ao«V 



2r 2 r 2 



= 



and 



- a u^ + 



n + (3k; + 



2 , 3/?o 



V0 = 0. 



(27) 



(2? 



Solving for n we have 



n = —(5 kl 



5ik r (3 
4r 



± 



c^fe 2 - 



16r 2 



+ 



2r 



1/2 



To find the real part of n, first write the bracketed term on the right of (|29|) as 

(a + frz) = (c + di) 1/2 , 



where 



and 



2,2 25Ay 2 /3 2 (k r Po\ ( 2 25^ 



16r 2 



<i = 3/c r aQ/2r 



h J y 1 16r 2 



(29) 



(30) 



(31) 



(32) 



2rh 2 ' 

and where (^) and (^) have been used. Note that a,b,c,d (|3TD, ( |32"D are all real. Solving 
for a and 6 we get 



1 



and 



a = ±^ r [c+{c 2 + d 2 ) 1 ' 2 ] 1 / 2 
v2 

^i^f-c+^ + rf 2 ) 1 / 2 ] 1 / 2 . 
v2 



(33) 
(34) 



Combining with (E3) we get for the real part of n: 



Re(n) ~ -ft** ± -L ^ 



2 25/i 2 / 4 /25/i r 
+ U 4 + 



1/2- 



16r 2 



16r 2 



1/2 



(35) 



where in calculating (c 2 +d 2 ), I have neglected terms of order q 2 k 2 /h 2 r 2 and q 4 /h 4 r 2 compared 
to a term of the order q A k 2 /h 4 , on the assumption that the relevant growth modes have 
k » 1/r (verified later) . 
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For q < hh/Ar, only the negative term in fl35|) survives and we have only decaying 
solutions. For q » 5h/4r we have 



Re(n) = ±k r p q/h — k 2 (3 . 



(36) 



Taking the postive sign, growth occurs for k r < q/h and the maximum growth rate occurs 
at k max = q/2h with maximum growth rate q 2 f3o/4h 2 . The growth is strongly dependent on 



A similar analysis shows that the imaginary part of n is given by 



Im(n) 



5k r p 1 / krpo 
M h , 



ir 



2hh 2 
16r 2 



25h 2 
16r 2 



1/2- 



1/2 



(37) 



for which only the first term on the right contributes when q >> 5h/4r. When q << 5h/4r 
the positive solution gives zero imagnary part, whilst the negative soultion gives Im(n) = 
-5k r f3 /2r. 



5. Discussion 

5.1. Planet formation 

The simple approach in the previous section turns out to provide vorticity growth with 
numbers comparable to what is required by planet formation. 

Studies of planet formation have invoked the idea of trapping dust particles in vortices, 
to provide the required agglomeration of material (Barge & Sommeria 1995; Tanga et al 1996; 
Hodgson & Brandenburg 1998 Godon & Livio 1999; Bracco et al. 1999 Chavanis 2000) In all 
of the above references, anti-cyclonic vortices (vortices opposed to the underlying rotation) 
are the only ones of interest for planet formation. Consider a vortex centered in the disc at 
radius r . Dust particles entering the vortex from r > r lose angular momentum and fall 
inward toward the vortex core, while those entering from r < r$ gain angular momentum 
and move outward, also toward the radius of the vortex core. While there is a limited size 
range of dust particles that couple to the vortex, the result is migration of material to the 
center for suitable anti-cyclonic vortices. 

To see that the discussion of the previous section leads to anti-cyclonic vortex growth, 
first note that ultimately, «o results from a vertical density gradient and the Coriolis force 
(hidden in the averaging rather than explcitly considered), so «o is positive in the northern 
hemisphere and negative in the southern hemisphere. Now the dominant growth term on 
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the right of fl22|) is the first term. That term is thus negative in the north and positive 
in the south. Thus the toroidal vorticity grows a negative contribution in the north and 
positive contribution in the south. This means that in (£23j), the toroidal velocity and thus 
the poloidal vorticity grow negative contributions (anti-cyclonic) in the north and also in 
the south. 

If we consider the case for which q ~ 1 » h/r, then we have k max ~ l/2h and n max = 
(3 /Ah 2 from the previous section. When the vortex reaches saturation, its characteristic 
evolution time would then be given by the oscillating part of the solution. From the q » 
5h/4r limit of ([3^]) we have 

2r 2r 

r v >2rh/f3 — — , (38) 

&ssCs il(y. ss \L 

where we have used k ~ k max . For a ss ~ 0.01 and h/r ~ 0.1, this limit would give > 300 
orbital periods. This is consistent with what is required to form planets (c.f. Godon & Livio 
1999). 

The saturated energy density associated with the mean field vortex should be no greater 
than that associated with the turbulent energy density. Using the fact that the growth time 
for the magneto-shearing instability driving the growth of turbulence in accretion discs is of 
order the dominant eddy turnover time, and is equal to the rotation period we have 
for the viscosity /3 Q = a ss c s h ~ v T /VL, where vt is the dominant turbulent speed. Thus 
vt ~ a ls c s and a vortex with the associated velocity for a mode of k r ~ q/h would have 

w ~ aH 2 c s k max ~ a l s i 2 n\q\ ~ O.lfi, (39) 

for q 2 » h 2 /r 2 , and a ss ~ 0.01. This is typical of that employed by by Godon & Livio 
(1999). 

Rewriting Godon & Livio's (1999) allowed range in size of surviving sub-sonic vortices 
we have 

- < A:" 1 < f^y (40) 

where the lower bound is a viscous length, the upper bound ensures that the vortex is not 
subject to destruction by shear, and we have used vt for the vortex speed. Using the relations 
for vt and (3$ above, we get the condition a]l 2 h < k~ l < al^frh) 1 / 2 . For k ~ h' 1 the lower 
inequality is always satisfied, and the upper inequality is satisfied for h/r < a 1 / 2 which is 
consistent with the range allowed at least in FU Orionis for which a ss could be even as large 
as 0.1 (Bell et al. 1995). 

These values are encouraging but it is also important to understand why a vortex core 
appears at a particular radius. For this simple-minded approach, the radial dependence of 
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the growth rate could enter through any hidden dependences in q. One could extremize 
(q — 5h/4r) to find the radius corresponding to the maximum. However, it is likely that 
the favorable radius of planet formation ultimately has much to do with the ambient dust 
concentration and size distribution as a function of radius in a stellar nebula (Barge & 
Sommeria 1995). 

It is very important to note that the overly simplified approach herein does not locate 
the vortices on azimuth, as the azimuth is averaged over. Thus for a given radius, only the 
net vorticity is measured. The location of the vortex, or the number of the vortices is not 
determined. 



5.2. Implications for accretion disc dynamos 

The growth of mean vorticity and mean magnetic field require the presence of a pseu- 
doscalar. It is interesting to compare the resulting pseudoscalars ao and a m by including 
the magnetic back-reaction to lowest order. 

Note that ao is proportional only to the negative of kinetic helicity whilst a m represents 
the residual between this and a current helicity contribution. (The relation between this and 
the form in Pouquet et al. (1976) is discussed in Field et al. 1999). This can in principle be 
tested in simulations, though one must be careful not to use periodic boundary conditions 
over the simulation box half-thickness over which the vertical averaging is done. In addition, 
it was noted in Brandenburg & Donner (1997) and Brandenburg (1999) that a rn seems to 
have the opposite sign than expected from just the kinetic helicity part. It was suggested 
that the shear may play the role of flipping the sign of a rising loop of field (see Brandenburg 
& Donner 1997). Here we can go a step further and note that it is actually the current 
helicity term that will be affected by that sign: imagine a seed toroidal field from which a 
loop rises and twists less than 1/4 turn, in the direction opposite to the underlying rotation. 
This would have positive current helicity. But now as the Keplerian shear takes over, the 
loop rotates past a 1/4 turn, and the sign of the current helicity changes. By contrast, the 
sign of the vorticity growth coefficient ao is the same for a Keplerian and non-shearing flow, 
since the kinetic helicity does not change sign. Thus in magneto-shearing simulations in 
which the turbulent magnetic energy slightly dominates the turbulent kinetic energy, the 
current helicity term may dominate a m . 
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5.3. Implications for accretion disc modeling 

The growth of vorticity as modeled herein may have implications for angular momentum 
transport in accretion discs. Indeed the presence of anti-cylconic vortices in the underlying 
flow transports angular momentum out of the ambient flow. Note however, that in the 
present treatment we have taken a base turbulent state that acts as a viscosity, so angular 
momentum is being transported also by the underlying turbulence. It is not entirely clear 
what the consequences of the vorticity growth would be for angular momentum transport, 
above and beyond their ability to agglomerate dust particles in the accretion flow. Such an 
additional effect of a needs to be calculated, since the pseudoscalar term is as naturally 
present as the viscosity term when integrating only over the half-thickness of the disc. 

Note that the generation of Rossby wave vortices has been considered when there is no 
turbulence present and leads to outward transport of angular momentum independently of 
any magnetic field. (Lovelace et al. 1999; Li et al 1999). 

Additional consequences of vortex generation for high energy accretion discs may include 
concentrating disc emission into collimated beams ( [Acosta-Pulido et al. 1990] , [Abramowicz 



1992[ ). X-ray iron line modeling could include the presence of such vortices to diagnose their 
presence. Also, Yoshizawa & Yokoi (1993) discussed how the interplay of magnetic field and 
vorticity can lead to generation of collimated large scale jets. 



6. Conclusions 

Standard axisymmetric Shakura-Sunyaev type turbulent accretion disc theory should 
always be regarded as mean field theory. When the vertical averaging is taken to be inte- 
gration over 1/2 the disc scale height, the presence of a pseudoscalar transport coefficient 
in addition to the usual scalar diffusion term should survive. This pseudoscalar term allows 
vorticity growth in each hemisphere even when mean quantities are only a function of radii 
there. The simplest growth rates, survival times and spatial scales seem to be consistent with 
that required by planet formation studies. The simple pseudoscalar term driving vorticity 
growth differs from that which drives mean magnetic field growth in that the latter, unlike 
the former, could have an opposite sign for cases of a rigid rotator vs. a sheared rotator such 
as an accretion disc. 

All of these points could be tested numerically, though one must be very careful about 
boundary conditions. Future analytic work can investigate the implications of including the 
pseudoscalar term in a generalization of the Shakura-Sunyaev/slim disc models. Azimuthal 
dependence must also be considered. 
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